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Abstract 

We study the dynamics of monodisperse hard elhpsoids via a new event-driven molecular dynam- 
ics algorithm as a function of volume fraction (j) and aspect ratio Xq. We evaluate the translational 
Dtrans and the rotational Drot diffusion coefficient and the associated isodiffusivity lines in the 
<j) — Xq plane. We observe a decoupling of the translational and rotational dynamics which gen- 
erates an almost perpendicular crossing of the Dtrans and Drot isodiffusivity lines. While the self 
intermediate scattering function exhibits stretched relaxation, i.e. glassy dynamics, only for large 
<j) and Xq ~ 1, the second order orientational correlator C2{t) shows stretching only for large and 
small Xo values. We discuss these findings in the context of a possible pre-nematic order driven 
glass transition. 

PACS numbers: 64.70.Pf,61.20.Ja,61.25.Em,61.20.Lc 



Particles interacting with only excluded volume interaction may exhibit a rich phase dia- 
gram, despite the absence of any attraction. Spherical objects, in equilibrium, present only 
a fluid and a crystal phase, while simple non-spherical hard-core particles can form either 



crystalline or liquid crystalline ordered phases [ll], as flrst shown analytically by Onsager 
for rod-like particles. Successive works have established detailed phase diagrams for sev- 
.„d-body ..apel. fl fl and Kave CanSed .Ke .ole of the e„.op. . t.e .an.tion 
between different phases. Less detailed information are available concerning dynamic prop- 
erties of hard-core bodies and their kinetically arrested states. In the case of the hard-sphere 
system, dynamics slows down signiflcantly on increasing packing fraction 0, and, when crys- 
tallization is avoided (mostly due to intrinsic sample polydispersity) , a dynamic arrested 
state (a glass) with extremely long life time can be generated. The slowing down of the 
dynamics is well described by mode coupling theory (MCT)i7]. On going from spheres to 
non-spherical particles, non-trivial phenomena arise, due to the interplay between transla- 
tional and rotational degrees of freedom. The slowing down of the dynamics can indeed 
appear either in both translational and rotational properties or in just one of the two. 

Hard ellipsoids (HE) of revolution [l,^] are one of the most prominent systems composed 
by hard body anisotropic particles. HE are characterized by the aspect ratio Xq = a/b 
(where a is the length of the revolution axis, b of the two others) and by the packing fraction 
(f) = TcXob^N/QV , where N is the number of particles and V the volume. The equilibrium 
phase diagram, evaluated numerically two decades ago j^, shows an isotropic fluid phase 
(I) and several ordered phases (plastic solid, solid, nematic N). The coexistence hues show 
a swallow-like dependence with a minimum at the spherical limit Xq = 1 and a maximum 
at Xq ^ 0.5 and Xq 2 (cf. Figure [I]). Application to HE(loj of the molecular MCT 
(MMCT) jlll . [l^ predicts also a swallow-like glass transition line. In addition, the theory 
suggests that for Xq ^ 0.5 and Xq ^ 2, the glass transition is driven by a precursor of 
nematic order, resulting in an orientational glass where the translational density fluctuations 
are quasi-ergodic, except for very small wave vectors q. Within MCT, dynamic slowing down 
associated to a glass transition is driven by the amplitude of the static correlations. Since 
the approach of the nematic transition line is accompanied by an increase of the nematic 
order correlation function at g = 0, the non-linear feedback mechanism of MCT results in a 



glass transition, already before macroscopic nematic order occurs 
rotational motions become hindered. 



10|. In the arrested state. 
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The recently evaluated random close packing line 13j for HE also exhibits a swallow- 



like shape. Although Kerr effect measurements in the isotropic phase of liquid crystals have 
given some evidence for the existence of two types of glass transitions [l4 ] (related to nematic 
phase formation and to cage effect respectively) almost nothing is known about the glassy 
dynamics of system forming liquid crystals in general and for HE in particular. 

We perform an extended study of the dynamics of monodisperse HE in a wide window of 



and Xq values, extending the range of Xq previously studied[15j. We specifically focus on 
establishing the trends leading to dynamic slowing down in both translations and rotations, 
by evaluating the loci of constant translational and rotational diffusion. These lines, in the 
limit of vanishing diffusivities, approach the glass-transition lines. We also study trans- 
lational and rotational correlation functions, to search for the onset of slowing down and 
stretching in the decay of the correlation. We perform event-driven (ED) molecular dynamics 
simulations, using a new algorithm [16|, which differently from previous algorithms [15|, [ITJ, 
relies on evaluations of distance between objects of arbitrary shape. We simulate a system 
of = 512 ellipsoids at various volumes V = in a. cubic box of edge L with periodic 
boundary conditions. We chose the geometric mean of the axis / = \^alP as unit of distance, 
the mass m of the particle as unit of mass (m = 1) and ksT = 1 (where ks is the Boltzmann 
constant and T is the temperature) and hence the corresponding unit of time is ^JmP /kBT. 
The inertia tensor is chosen as = ly = 2mr^/5, where r = min{a, b}. The value of the 
component is irrelevant [18!], since the angular velocity along the symmetry (z-) axis of the 
HE is conserved. We simulate a grid of more than 500 state points at different Xq and as 
shown in Fig. [H To create the starting configuration at a desired 0, we generate a random 
distribution of ellipsoids at very low and then we progressively decrease L up to the de- 
sired 0. We then equilibrate the configuration by propagating the trajectory for times such 
that both angular and translational correlation functions have decayed to zero. Finally, we 
perform a production run at least 30 times longer than the time needed to equilibrate. For 
the points close to the I-N transition we check the nematic order by evaluating the largest 
eigenvalue S of the order tensor Q jl^, whose components are: 

Qaf3 = J]((Ui)a(Ui)^) - -5a^f3 (l) 

i 

where a/3 G {x, y, z}, and the unit vector {ui(t))a is the component a of the orientation (i.e. 
the symmetry axis) of ellipsoid i at time t. The largest eigenvalue S is non-zero if the system 



is nematic and if it is isotropic. In the following, we choose the value S = 0.3 as criteria 
to separate isotropic from nematic states. From the grid of simulated state points we build 
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FIG. 1: Grid of state points simulated (crosses and red-filled circles) and relevant boundary lines of 



I_ first order phase transitions 
. Solid lines are coexistence 



coexistence regions. Long-dashed curves are coexistence curves of a 
in the phase diagram of HE evaluated by Frenkel and Mulder (FM) 
curves for the I-N transition of oblate and prolate ellipsoids, obtained analytically by Tijpto-Margo 
and Evans (TME). Dotted lines (TME-2) are coexistence curves of prolate ellipsoids for the I-N 
transition, taken from j^. 

a corresponding grid of translational (Dtrans) and diffusional (Drot) coefficients, defined as: 
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where A$j = cuidt, Xj is position of the center of mass and tOi is the angular velocity 
of ellipsoid i. By proper interpolation, we evaluate the isodiffusivity lines, shown in Fig. 
[21 Results show a striking decoupling of the translational and rotational dynamics. While 
the translational isodiffusivity lines mimic the swallow-like shape of the coexistence between 
the isotropic liquid and the crystalline phases (as well as the MMCT prediction for the 
glass transition [10|]), rotational isodiffusivity lines reproduce qualitatively the shape of the 
I-N coexistence. As a consequence of the the swallow-like shape, at large fixed 0, Dtrans 
increases by increasing the particle's anisotropy, reaching its maximum at Xq ^ 0.5 and 



Xo ~ 2. Further increase of the anisotropy results in a decrease of Dtrans- For all Xq, an 
increase of (f) at constant Xq leads to a significant suppression of Dtrans, demonstrating that 
Dtrans IS coutroUed by packing. The iso-rotational lines are instead mostly controlled by 
Xq, showing a progressive slowing down of the rotational dynamics independently from the 
translational behavior. This suggests that on moving along a path of constant Dtrans, it is 
possible to progressively decrease the rotational dynamics, up to the point where rotational 
diffusion arrest and all rotational motions become hindered. Unfortunately, in the case of 
monodisperse HE, a nematic transition intervenes well before this point is reached. It is 
thus stimulating to think about the possibility of designing a system of hard particles in 
which the nematic transition is inhibited by a proper choice of the disorder in the particle's 
shape/elongations. We note that the slowing down of the rotational dynamics is consistent 
with MMCT predictions of a nematic glass for large Xq He[i^, in which orientational 
degrees of freedom start to freeze approaching the isotropic- nematic transition line, while 
translational degrees of freedom mostly remain ergodic. To support the possibility that the 
slowing down of the dynamics on approaching the nematic phase originates from a close-by 
glass transition, we evaluate the self part of the intermediate scattering function Fseif 



and the second order orientational correlation function C2(t) defined as [l5| C2(t) = 
{P2{cos 6 (t))) , where P2{x) = (3x^ — l)/2 and 6{t) is the angle between the symmetry 
axis at time t and at time 0. The C2{t) rotational isochrones are found to be very similar 
to rotational isodiffusivity lines. 

These two correlation functions never show a clear two-step relaxation decay in the entire 
studied region, even where the isotropic phase is metastable, since the system can not be 
significantly over-compressed. As for the well known hard-sphere case, the amount of over- 
compressing achievable in a monodisperse system is rather limited. This notwithstanding, a 
comparison of the rotational and translational correlation functions reveals that the onset of 
dynamic slowing down and glassy dynamics can be detected by the appearance of stretching. 
Fig. [3] contrasts the shape of F^eif, evaluated at g = qmax, where qmax is the q corresponding 
to the first maximum of the center-of-mass static structure factor, and C2{t) at = 0.50 
for different Xq values with best-fit based on an exponential (~ exp[— t/r]) and a stretched 
exponential (~ exp[— (t/r)^]) decay. As a criteria to avoid including in the fit the short-time 




(4) 
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FIG. 2: Isodiffusivity lines. Solid lines are isodiffusivity lines from translational diffusion coefficients 
D trans and dashed lines are isodiffusivities lines from rotational diffusion coefficients Drot- Arrows 
indicate decreasing diffusivities. Left and right arrows refer to rotational diffusion coefficients. 
Diffusivities along left arrow are: 1.5, 0.75, 0.45, 0.3, 0.15. Diffusivities along right arrow are: 1.5, 
0.75, 0.45, 0.3, 0.15, 0.075, 0.045. Central arrow refers to translational diffusion coefficients, whose 
values are: 0.5, 0.3, 0.2, 0.1, 0.04, 0.02. Thick long-dashed lines are FM and TME coexistence lines 
from Fig(T] 



ballistic contribution, we limit the time- window to times larger than t*, defined for -Fsez/ and 
6*2 as the time at which the autocorrelation function of center-of-mass velocity v (^^(t) = 
14 Z]i(vi(i)vi(0)) ) and of angular velocity respectively (0<^<^(t) = Z]j(^i(^)^i(0))) reaches 
1/e of its initial value. We note that -Fsei/ shows an exponential behaviour close to the I-N 
transition (Xq = 3.2, 0.3448) on the prolate and oblate side, in agreement with the fact 
that translational isodiffusivities lines do not exhibit any peculiar behaviour close to the I-N 
line. Only when Xq ~ 1, Fseif develops a small stretching, consistent with the minimum 
of the swallow-like curve observed in the fluid-crystal line 201. 2l[ |. in the jamming locus 



as well as in the predicted behavior of the glass line for HE lOj and for small elongation 



dumbbells 
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231]. Opposite behavior is seen for the case of the orientational correlators. 



C2 shows stretching at large anisotropy, i.e. at small and large Xq values, but decays within 
the microscopic time for almost spherical particles. In this quasi-spherical limit, the decay 
is well represented by the decay of a free rotator[24]- Previous studies of the rotational 
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FIG. 3: Shape of Fgeif and C2 at (p = 0.50 for different Xq. Symbols are data from MD simulations. 
Solid lines are fits to exponential functions, while long-dashed lines are fits to stretched exponentials 
(/? is the stretching parameter) . t* is the time at which correlation functions (p^ and ; for Fgei / 
and C2 respectively, reach 1/e of their initial values. Top: Prolate ellipsoids with Xq = 3.2, 
C2 shows a significant stretching while F^^if decays exponentially. Center: Xq = 1.0 for Kje// 
and Xq = 1.1 for C2, the dashed line is the theoretical decay of a free rotator [Cl^t) = 
1 - ^^exp[-t'^/Tj]^{t/Tf), where rj = l/<pujuj{0) and = exp[x'^]dx) . Bottom: Oblate 
ellipsoids with Xq = 0.348. 



dynamics of HE [151] did not report stretching in C2, probably due to the smaller values of 
Xq previously investigated and to the present increased statistic which allows us to follow 
the full decay of the correlation functions. 

Fig. [3] clearly shows that C2 becomes stretched approaching the I-N transition while Fseif 
remains exponential on approaching the transition. To quantify the amount of stretching 
in C2 we show in Fig. H] the Xq dependence of r and (3 for three different values of 0. In 
all cases, slowing down of the characteristic time and stretching increases progressively on 
approaching the I-N transition. It is interesting to observe that the amount of stretching 
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appears to be more pronounced in the case of oblate HE compared to prolate ones. A similar 
(slight) asymmetry between oblate and prolate HE can be observed in the lines reported in 
Figure 2. 

In summary, we have shown that clear precursors of dynamic slowing down and stretching 
can be observed in the region of the phase diagram where a (meta) stable isotropic phase 
can be studied. Despite the monodisperse character of the present system prevents the 
possibility of observing a clear glassy dynamics, our data suggest that a slowing down in the 
orientation degree of freedom — driven by the elongation of the particles — is in action. The 
main effect of this shape- dependent slowing down is a decoupling of the translational and 
rotational dynamics which generates an almost perpendicular crossing of the Dt^ans and D^-ot 
isodiffusivity lines. This behavior is in accordance with MMCT predictions, suggesting two 
glass transition mechanisms, related respectively to cage effect (active for 0.5 ^ Xq ^ 2) and 
to pre-nematic order (Xq ^ 0.5, Xq ^ 2) [lO|. It remains to be answered if it is possible to 
find a suitable model, for example polydisperse in size and elongation, for which nematization 
can be sufficiently destabilized, in analogy to the destabilization of crystallization induced 
by polydispersity in hard-spheres. We acknowledge support from MIUR-PRIN. We also 
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FIG. 4: [3c2 and obtained from fits of C2 to a stretched exponential for (j) = 0.40, 0.45 and 

0.50. Top: as a function of Xq. Bottom: [3c2 as a function of Xq. The time window used for 
the fits is chosen in such a way to exclude the microscopic short times ballistic relaxation (see text 
for details). For 0.588 < Xq < 1.7 the orientational relaxation is exponential. 
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